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1. Introduction 

In a series of papers®^ we have argued that the effective cosmological constant 
A of the present universe, interpreted as a remnant vacuum energy density, may 
be due to the imprint of ultraheavy particles on the Hubble expansion of the early 
universe. The discussion is in the framework of the so-called g-theory which describes 
the evolution of the macroscopic gravitating vacuum energy density pv [ 9 ] due to a 
microscopic conserved relativistic vacuum variable q (the original papers are Refs. SI 
[S]and a one-page review appears as App. A in Ref. IH]). In view of the upcoming Run 
2 of CERN’s Large Hadron Collider (LHC), it may be timely to review the main 
assumptions of the argument and to clarify the meaning of the predictions. Natural 
units will be used with h = 1 and c = I. 

At this point, it may already be worthwhile to present the basic equation for 
the remnant vacuum energy density in a flat Friedmann-Robertson-Walker (FRW) 
universe with cosmic time t, 

A = \im pv{t) =rvoo /{Epf , (I.l) 

t—¥(yD 

where rvco is a nonnegative number, M the mass scale of the hypothetical new 
particles, and Ep the reduced Planck energy, 

Ep = y/l/(87rGjv) R:! 2.44 X 10^® GeV. (1.2) 
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Inverting dm gives the following expression for the new mass scale M: 

/ in-3 \ 1/® / A 1/4 \ 

M = A'/* 5.56 Tev(_) { ^;^) , (1-3) 

where the numerical value used for A follows from Table 2.1 in Ref. [71 

Taking Eq. (EB as it stands, the number rvoo can be interpreted as an ef¬ 
ficiency factor for producing a remnant vacuum energy density given the energy 
scales involved, M and Ep. Remark that the parametric dependence of EB has 
already been discussed by Arkani-Hamed et al.,^but without a convincing theory 
of how the vacuum energy density evolves (it is here that g-theorjJ^I^ is supposed 
to take over). Still, the authors of Ref. [8] discuss persuasively the role of (11.11) for 
the so-called triple cosmic coincidence puzzle: why are the orders of magnitude of 
the energy densities of vacuum, matter, and radiation approximately the same in 
the present Universe? 

Taking the point of view that the effective cosmological constant A of the present 
universe has been measured,!^ Eq. (11.31) can be read as a prediction of the mass scale 
M, provided the “efficiency factor” rvoo is known. 

The task, then, is to calculate the pure number rvca entering (II. ip . This cal¬ 
culation is extremely difficult and, up till now, there is only an approximate phe¬ 
nomenological description availableIn the present paper, we try to simplify the 
discussion as much as possible, in order to highlight the crucial assumptions of 
the argument. For definiteness, we assume q to come from the field strength of a 
three-form gange field,® so that q has mass dimension 2. 

2. Setup 

2.1. K-freezing model 

In the framework of g-theory,®® the analysis of Ref. [T] has shown that the sudden 
presence of ultraheavy particles (possibly created by a phase transition with decay 
afterwards) perturbs the Hubble expansion and kicks the vacuum energy density 
pv{t) away from zero to a small positive value. This process occurs at a cosmic age 
of order 

4ick = Ep/M^ = M-i, (2.1) 

where the last expression has been written in terms of the energy-density hierarchy 
parameter 

^={Ep/M)\ (2.2) 

In the same way as the authors of Ref. [HI we initially do not worry about what 
physics stabilizes the large hierarchy between M ^ TeV and Ep ^ 10^^ TeV. In 
other words, we leave aside the well-known hierarchy problem and simply try to 
determine the mass scale M from cosmology, for the given value of Ep. 
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The fundamental issue is the freezing of the vacuum energy density created by 
the kick. As realized in Ref. 121 this freezing can be modeled by a time-dependent (or, 
better, temperature-dependent) gravitational coupling; see also Sec. 12.21 for further 
comments. Specifically, we take a Brans-Dicke-type term in the action density, 

C^rav = K[q,^R[g], (2.3) 

where $ stands for one or more of the new matter fields. Recall that the standard 
Einstein theory has an action-density term Kq R[g] with constant Kq = l/(167rGjv)- 
In a flat FRW universe with cosmic time t, we now assume the following simplified 
behavior: 


K[q, t] = q{t)/2 -k (t) [qo/2 - q{t)/2] , (2.4a) 


dK{t) = 9{t — Ik) , 


in terms of the standard stepfunction 


e{t) 


1 for t > 0 , 
0 for t < 0 . 


(2.4b) 


(2.5) 


Observe that q^ in (I2.4al) corresponds to the constant equilibrium value of the vac¬ 
uum variable q and that the inverse of qo gives Newton’s constant, Gjv = 1 /(Stt^q). 

We also let the coupling constant of the ultraheavy matter component be con¬ 
trolled by another stepfunction. 


g^(t) = g^ 9g{t), 


(2.6a) 


9g{t) = 9{t-tg). 


(2.6b) 


This time-dependent coupling constant generates an ultraheavy matter component, 
even if it is not present initially (t < tg). At a later moment, the ultraheavy particles 
decay, which can be modelled by letting g'^{t) drop to zero again and by having a 
constant decay constant instead. For the kick mechanism to operate (that is, ry 
first kicked away from zero and then frozen), the following inequality is required: 


tg <tK ■ 


(2.7) 


Moreover, the kick mechanism only makes sense if these two timescales are of the 
same order of magnitude (see Sec. H] for further discussion). 

Next, define dimensionless variable^ by use of the energy scales M and Ap, 
together with the auxiliary hierarchy parameter ^ from (12.21) . The dimensionless 
cosmic time, in particular, is defined by t = t/tkick with tkick from (j2.1ll . The rele¬ 
vant dynamic variables are the dimensionless Hubble parameter /i(t), the rescaled 
relative q-parameter shift x{t), the rescaled dimensionless energy density rMiix) 
of the ultraheavy particles (called type-1), and the rescaled dimensionless energy 
density rM 2 {T) of the massless particles (called type 2). The type-1 particles are, for 
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definiteness, assumed to be bosons. In a finite temperature context without chem¬ 
ical potential, the type-1 bosons have an equation-of-state parameter wmi, which, 
depends only on M and T. The massless type-2 particles have the equation-of-state 
parameter wm 2 = 1/3 and the effective number of degrees of freedom Nes ,2 = 10^, 
in order to represent the lighter particles of the Standard Model (see Sec. V.B of 
Ref. m for further discussion). With the temperature T obtained from the energy 
density Pm 2 = 2 (’’’/SO) T"*, it is then possible to write wmi in terms of M and 

Pm 2 and the resulting expression is denoted wmi', see Sec. A2 of Ref. [2] for details. 
Finally, we define the combination 

kmi = 1 — 3 Wmi , (2.8) 

which has been found to drive the kick of the vacuum energy density.® 

The ordinary differential equations (ODEs) for the four dynamic variables ar^ 


(1- 

Ok) 

3h + Qh'^ — X 

+ 0k 

6hh - Tmi - rMa] = 0, 

(2.9a) 

rMi 

+ (4 

— KMi) hr Ml 

= +5^ ^M2 — g^ tmi — tmi, 

(2.9b) 

f’M2 

+ 4h 

rM2 = -g^ rM2 + g"^ 

Tmi + A^ Tmi, 

(2.9c) 

(1- 

Ok) 

Shx/^ + 

x/^- 

^x^/i^i) + Tmi + rM2 - 3 


+0k 

i = 

0, 



(2.9d) 


with 9k given by (|2.4bl) and by (12.61) . both in term of the dimensionless cosmic 
time T (the overdot in these ODEs denotes differentiation with respect to r). The 
source terms in Eqs. (I2.9bl) and (I2.9cl) are somewhat different compared to those 
of (3.3b) and (3.3c) in Ref. [3] Equations (I2.9al) and (I2.9dl) correspond to (3.3a) 
and (3.3d) of Ref. |3] but are slightly rewritten (actually, these rewritten ODEs 
were already used for the numerics presented in Ref. [3]). Note that (I2.9al) in the 
9k = ^ phase (late times) gives the derivative of the standard Hubble equation 
— tmi — »'M 2 ~ 'k’v/^ = 0 with a term ry = (1/2) = const., according to 

(I2.9dp for 9k = 1- 

An exact solution of the ODEs (12.9p for = 0 is given bj® 


h(r) = 1/(2 r), 

(2.10a) 

x { t ) =0, 

(2.10b) 

TMiir) = 0, 

(2.10c) 

t’M 2 (t) = 3 [/i(t)]^ = rM 2 ('ro) [a(ro)/a(r)]'‘, 

(2.10d) 


where the last expression uses the scale factor a(r), in terms of which /i(t) is defined 
by a(r)/a(T). The exact solution (I2.10|) holds for both phases, the early one with 
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Ok = 0 and the late one with Ok = 1- Incidentally, the dimensionless cosmic age 
T = 0.27 corresponds to having T M. 

The goal, now, is to study the kick of ry (r) by the sudden presence of ultraheavy 
particles, which, in our model, results from a nonzero coupling constant g^ir) for 
T > Tg > Tijiin. For this purpose, the ODEs with ^^(t) from (12.61) are to be solved 
using the following boundary conditions: 


^('Tmin) 

= 1/(2 Tmin) , 

(2.11a) 

^('^min) 

= 0, 

(2.11b) 

7'Ml('rmin) 

= 0, 

(2.11c) 

?'M2('rmin) 

= 3 [/l(rmin)]^ , 

(2.lid) 


which matches the solution (I2.10|) of the earlier phase r < 

For later use, we can mention that the ^ = oo equations have been discussed in 
Ref. 13] and that an analytic result has been obtained for the vacuum energy density 
without iti-freezing effects. 


ry(r) 


(^=oo , no i<’-freezing) 


KMi(r) rMi(r) 


1 2 


( 2 . 12 ) 


But (]2.12p still needs two ^ = oo ODEs to be solved, in order to obtain the explicit 
functions Umi (t) and tmi (t) . 


2.2. Dissipation model 

In this article, we primarily model the freezing of ryir) by taking a particular 
time-dependent K[q,t] and by remaining entirely within g-theory which describes 
reversible processes. But, as argued in Sec. IV of Ref. (TJ the freezing of rv(r) 
may very well be due to quantum-dissipative effects, that is, irreversible processes.^ 
Concretely, very low-energy gravitons (and possibly neutrinos) may be responsible 
for the dissipation.ESl Specifically, the energies of these particles must be of the 
order of meV, which value traces back to (EB. In any case, a phenomenological 
description of this quantum dissipation has been given in Ref. |TJ which we will now 
briefly review. 

The following model equatiorP can be used: 


= -r(T) v(r) - rv.o(r) 


(2.13) 


which recalls the standard description of bulk-viscosity effects in fluid mechanics 
[see, in particular, Eq. (78.1) of Ref.|9]. The interpretation of (I2.13|) is that rv,o(T) is 
the “bare” vacuum energy density driven by the kick from the ultraheavy particles 
and that r(r) > 0 is the rate at which the “excess” vacuum energy density is 
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dissipated into particles. Equation (12.131) has an exact solution^ 

= J dr'r(T')rv,o(T') exp -j dr" T{t") , (2.14) 

for boundary condition py(0) = 0, which holds for times r well after the Planck 
era. 

In the next section, we will obtain a numerical estimate of the asymptotic value 
of based on (12.131) . by use of the following approximations: 


If- 

»'V,o(t) = - KMliT)rMl{T) ■ 


(2.15a) 


r(r) = 7/r2 , (2.15b) 

where the first approximation becomes exact in the limit ^ > oo, according to (j 2 . 12 |) . 
and the second approximation is purely illustrative. Taking a relatively large value 
for 7 forces r^'®®(r) to follow rv,o('7‘), according to (12.1311 or (12.141) . A small enough 
finite value of 7 allows ry®®(T) to reach a nonzero asymptotic value. 

Remark that (12.131) is only part of the whole story, as it does not specify the 
detailed changes in the matter energy densities PM,n{T) corresponding to the change 
of But this partial description is perfectly valid during the kick phase, as the 

amount of energy carried by the vacuum then is negligible compared to that of the 
ponderable matter by a factor of order ^ ~ 10®^, according to (12.211 for M ~ 10 TeV 
(see also Fig. 1 of Ref. IH]). 


3. Numerical results 

The numerical solutions of the ODEs (j2.9ll with boundary conditions (12.111) are 
readily obtained. We present numerical solutions in Figs. [T] and [2] for two choices 
of coupling constants ( 5 , A), which give a ratio rMi( 0 . 3 )/rM 2 ( 0 . 3 ) of order 1/100 
and 1, respectively. In an equilibrium context, these energy-density ratios would 
translate into a degrees-of-freedom ratio 1 /A"eff, 2 = 1/100 for case 1 and a ratio 
Aeff, i/NeS ,2 = 100/100 for case 2. Note that the small oscillations on ry = (1/2) 
in the top-right panels of the figures disappear for even larger values of the hierarchy 
parameter, ^ ^ 10^. At this moment, it may be useful to recall the analytic ^ = 00 
result (I2.12|) . This analytic result suggest to consider the quantity 1/8 (kmi tmi) 
shown in the third bottom-row panels of the figures, which is indeed smoother than 
the numerical result ry{T) shown in the top-right panels. 

The ODEs (12.91) result from the AT-freezing model and the main result of the 
numerical calculation is the asymptotic plateau of the vacuum energy density ry (t) 
as shown in the top-right panels of Figs.[I]and[2] For comparison, we also show, in the 
bottom-right panels, the vacuum energy density resulting from the dissipation model 
(j2.13l) with approximations p.l5|l and an appropriately chosen 7 value. Similar 
results are obtained with r(r) Ansdtze that drop to zero faster than 1 /r^. An 
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Figure 1. Numerical solution of the dimensionless ODEs II2.9II with an equation-of-state function 
= 1 —3 uJmi a.s defined in Sec. A2 of Ref. [2] The panels are organized as follows: the four 
basic dynamic variables ^(t), r^ii'^)^ x(t) are shown on the top row and secondary 

or derived quantities on the bottom row. The dashes lines in certain panels refer to the second 
quantity listed in the respective panel label, for example, the dashed line in the bottom-left panel 
corresponds to 9k • The main result is the nonzero remnant value of the dimensionless gravitating 
vacuum energy density ry = x^/2 shown in the top-right panel. The bottom-right panel shows, for 
comparison, the dimensionless vacuum energy density from quantum-dissipative effects, as modeled 
by II2.13II with approximations (12.1511 . The model parameters are Neff, 2 ^ 5^? '^g-, '^K-, 7} = 

{10^, 10^, 1, 10^, 0.2, 0.4, 1/5} and the ultraheavy type-1 particles are assumed to be bosons 
(similar results are obtained for type-1 fermions, with somewhat lower values for ry and by 
approximately 20%). The ODEs are solved over the interval [rmixn "Tnax] = [0.02, 0.8] with the 
following boundary conditions at r = Tmin = 0.02: {x, h, a, ^M 2 } = {0, 25, 1/10, 0, 1875}. 
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Figure 2. Same as Fig. ^ but with model parameters A^} = {64, 1}. These coupling con¬ 
stants make for a ratio rM\l'^M 2 1 at r = {jg tk)I‘^ = 0.3, whereas the coupling constant of 
Fig. ^ give a ratio of order 1/100. 


exponential tail, for example, has been used in a previous version of the present 
paper [arXiv:1503.03858vl]. 

For the case of Fig. [T] (‘Weff, i = 1 ”) with an ry peak of order 2 x 10“^, the 
frozen asymptotic value ry^ obeys the upper bound 


(case-1) 


< max [r'y(r)] 


(case-1) 


10 


-5 


QO 


( 3 . 1 ) 
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which translates into the following lower-bound on M from IE3: 

(case-1) 


M 


> 10 TeV. 


(3.2) 


For the case of Fig. [5] (‘Weg^i = 10^”) with an ry peak of order 7 x 10 the 
asymptotic value rvoo obeys the upper bound 


rvc 


(case-2) 


< max [rv(T)] 


(case-2) 


10 


-1 


which gives 


M 


(case-2) 


> 3 TeV. 


(3.3) 


(3.4) 


4. Conclusion 

The results (13.21) and (13.41) from the Til-freezing model of Sec. I2.1l set the mass scale 
of the hypothetical new particles. If the underlying physics is able to relate the 
time scale tg of ultraheavy (type-1) particle creation and the time scale Ik of the 
change in the effective gravitational coupling, these inequalities could be replaced 
by rough equalities and be all the more convincing. Perhaps such a single physical 
process is similar to the one of slow-roll particle production discussed in the context 
of inflation models.!^ In our case, the slow-role phase must be relatively short. More 
importantly, such a process must not reinstate the cosmological constant problem 
(with a new scale of order M^) and must, therefore, include the evolution of q. 
Particle creatiorP^ is also crucial for quantum-dissipative effects of the vacuum 
energy density, as illustrated by the alternative model of Sec. 12.21 

As (13.21) and (13.41) are only lower bounds, we cannot predict that the mass scale 
M must necessarily be in reach of the LHC with a 13 TeV center-of-mass energy 
(for a possible high-luminosity upgrade, see Ref. (E]). Perhaps a next-generation 
proton-proton collider with 50 — 100 TeV center-of-mass energjff^E^ jg needed. For 
the moment, we can only adopt a “wait-and-see” attitude 0 If multi-TeV particles 
are discovered at the LHC or a next-generation collider (and, admittedly, this is a 
big ‘if’), an exciting prospect may be that they provide a new window on the early 
universe and fundamental physics. 
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^As to the properties of the hypothetical new particles, we remain entirely agnostic. One possibility 
involves ultraheavy composite scalars from top-condensation models for the light Higgs,but 
many other models exist. 
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